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We consider the effect of non-standard neutrino interactions (NSI, for short) on the propagation of 
neutrinos through the supernova (SN) envelope within a three-neutrino framework and taking into 
account the presence of a neutrino background. We find that for given NSI parameters, with strength 
generically denoted by Sij, neutrino evolution exhibits a significant time dependence. For {ettI 
10""^ the neutrino survival probability may become sensitive to the 'O23 octant and the sign of e^r. 
In particular, if Ett ^ 10^^ an internal /-resonance may arise independently of the matter density. 
For typical values found in SN simulations this takes place in the same dense-neutrino region above 
the neutrinosphere where collective effects occur, in particular during the synchronization regime. 
This resonance may lead to an exchange of the neutrino fluxes entering the bipolar regime. The 
main consequences are (i) bipolar conversion taking place for normal neutrino mass hierarchy and 
(ii) a transformation of the flux of low-energy p^, instead of the usual spectral swap. 

PACS numbers: 13.15.-|-g, 14.60.Lm, 14.60.Pq, 14.60.St, 97.60.Bw 



I. INTRODUCTION 

Current neutrino data imply that neutrino have mass. 
Indeed reactor data by the KamLAND collaboration 1| 
and data from accelerator neutrino experiments K2K [2| 
and MINOS Q not only confirm the neutrino flavor con- 
version discovered in the study of solar 0-0] and atmo- 
spheric neutrinos @, Q but also indicate that the under- 
lying neutrino flavor conversion mechanism is oscillatory 
in both cases [l3|- An updated review of the current 
status of neutrino oscillations is given in Refs. |11|. 

Theories of neutrino mass [l^ . [l3| typically require that 
neutrinos have non-standard four-Fermi interactions ^ 
which, for short, we denote as NSI [l5-17|. These are 
natural outcome of neutrino mass models and can be of 
two types: flavor-changing (FC) and non- universal (NU). 
For example, generic seesaw-type models lead to a non- 
trivial structure of the lepton mixing matrix character- 
izing the charged and neutral current weak interactions 



implying an effective non-unitary form for the mixing ma- 
trix describing neutrino oscillations [l^ . However, the 
expected magnitude of the unitarity violation and of the 
four-Fermi NSI effects is rather mo del- dependent. In the 
simplest high-scale seesaw models these are all negligi- 
ble. In contrast, the 4-Fermi NSI effects induced by the 
charged and neutral current gauge interactions may be 
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sizeable in low-scale seesaw schemes 

Alternatively, non-standard neutrino interactions may 
arise in models where neutrino masses are radiatively in- 
duced by low-scale loop effects 26l428l|. or directly "cal- 
culable" by renormalization group evolution [2^ in some 
supergravity models. 

It is important to realize that the strengths of non- 
standard interactions need not be suppressed by the 
smallness of neutrino masses. Indeed, they may be very 
relevant for the propagation of supernova (SN) neutri- 
nos and give rise to a novel type of resonant neutrino 
flavor conversion mechanism that can take place even in 

3^ 2 Therefore we 



the limit of massless neutrinos 
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Other non-standard neutrino properties such as electromagnetic 
transition moments may also be present [14| |. 



^ Similarly flavor and CP-violating effects can survive in the limit 
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argue that, in addition to the precision determination 
of the oscillation parameters, it is necessary to test for 
sub-leading non-oscillation effects that could arise from 
non-standard neutrino interactions in upcoming neutrino 
experiments [3l|, |32| . Especially because NSI effects can 



33 



3- 



in some cases fake genuine mixing effects 

Here we concentrate on the impact of non-standard 
neutrino interactions on SN physics. The main motiva- 
tion of the work is to reexamine the effect of non-standard 
neutrino interactions on supernova when the neutrino 
self- interaction is taken into account. The extreme con- 
ditions under which neutrinos propagate, since they are 
created in the SN core until they reach the Earth, may 
lead to strong matter effects. It is known that, in partic- 
ular, the effect of small values of the NSI parameters can 
be dramatically enhanced in the inner strongly dclcp- 
tonized regions 17, 30, 35, [sCL On the other hand, it 
has been recently pointed out 37H6lj that in this region 
the large neutrino background could itself affect neutrino 
propagation significantly. In this paper we analyze the in- 
terplay between these two in principle coexisting effects, 
the NSI and the neutrino self-interaction, and the result- 
ing consequences for neutrino evolution through the SN 
envelope. 

The paper is organized as follows. In Sec. [Tl] we de- 
scribe the equations of motion that govern neutrino prop- 
agation taking into account both NSI and self interac- 
tion effects. In Sec. IIIII we study the interplay between 
these two effects and analyze the conditions for different 
regimes of evolution to take place. In Sec. IIVI and |Vl 
we describe in detail the salient features characterizing 
neutrino evolution in the different regions defined by the 
conditions given in the previous section. Finally, before 
summarizing in Sec. IVIIl we give, in Sec. IVIl a time de- 
pendent discussion of the studied effects. 



II. EQUATIONS OF MOTION 

The equations of motion (EOMs) of the neutrinos trav- 
eling through the SN envelope can be written as 



(1) 



where pp and pp represent the density matrices describ- 
ing each (anti)neutrino mode. The diagonal entries are 
the usual occupation numbers whereas the off-diagonal 
terms encode phase information. The Hamiltonian for 
neutrinos has the form 

d3q 



(2) 



Hp = r!p+V+V2GF j (f?q - ^q) (1-Vq-Vp 

For antineutrinos the only difference is fip — > — ilp. 

The first term stands for the matrix of vacuum os- 
cillation frequencies, flp = C/diag(m^, TO2: "^D/^IpI^^^ 
in the weak basis, where U is the three-neutrino lep- 
ton mixing matrix 12| in the unitary approximation 
and PDG convention [62[, with no CP phases. We 
use Am^i = m| - = 7.65 x 10"^ eV^, |Am§i| = 
|m§ - ml\ = 2.40 x lO"^ gy^, sin^ i?i2 = 0.3, as ob- 



tained in e.g. Ref [ll|. We consider also sin^ 1^13 = 10 



and three different values for '62^ in the allowed range 
at 3(7, sin^ 1^23 = 0.4, 0.5 and 0.6, because our results 
depend sensitively on z?23- Given the values of Am^ one 
can define the two associated vacuum oscillation frequen- 
cies: wh = /S.m\i/2E and wl = ^'m\i/2E, which, in 
the case of neutrinos with = 20 MeV, take the values 
wh = 0.3 km^^ and wl ~ 0.01 km^^. In the top panel 
of Fig. [T]we represent wjj and ljl for energies typical in 
SNe, between 5 MeV and 50 MeV, as yellow and light 
blue horizontal bands, respectively. 

The second term of the Hamiltonian accounts for the 
interaction of neutrinos with matter and can be split into 
two pieces. 



V = V, 



std 



(3) 



The first term, Vgtd, describes the standard interaction 
with matter and can be represented in the weak basis 
by Vstd = V2GYnBd\&g{Ye,Q,Yf), with Gp the Fermi 
constant, ns the baryon density, and Ye = nf^/ns stands 
for the electron fraction. We consider the following stan- 
dard weak potential 



v,td = A(r)diag(ye,o,y;^), 



with 



A(r) = Ao 



R 



(4) 



(5) 



of masslcss neutrinos Il9l [20II . 



In the following we assume i? = 10 km. In the top panel 
of Fig.[T]wc show two A(r) profiles for Ai^o = 5x 10^ km~^ 
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FIG. 1: Top panel: A(r), profiles as defined in Eq. ((5)1, for 
Ai,o = 5 X 10^ km-i (Ai) and A2,o = 4 x 10*' km"^ (A2) 
in solid red and blue lines, respectively; A2(r), as given by 
Eq. pi|l . is shown in blue dotted lines; pi(r), introduced in 
Eq. (|13p . for /lo — 7 x 10^ km~^ in black dashed lines; the 
vacuum oscillation frequencies lju (yellow band), ujl (light 
blue band), and cjJJ", defined in Eq. (|17|l . for Ett ~ and 0.1 
(green bands), for energies between 5 and 50 MeV. The posi- 
tion of the synchronization and bipolar radii are also shown. 
Bottom: radial dependence of Ye, given by Eq. ((Gjl, for two set 
of parameters: a = 0.24, b = 0.165, ro = 50 (1.2 x 10^) km, 
and = 5 (3x 10^) km, for Y^" (Y^). The horizontal magenta 
band represents the YJ leading to an internal /-resonance for 

Err < 0.1. 



and A2,o = 4 x 10^ km ^ denoted by Ai and A2, corre- 
sponding to typical early and late time profiles, respec- 
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tively 

The first element in V'"^ represents the charged cur- 
rent potential and is proportional to the electron fraction. 
Ye- According to the SN models, Ye is characterized by 
a transition from a few % in the inner most deleptonized 
layers until values around 0.5 in the outer envelope. Fol- 
lowing Ref. [36[ we parametrize it phenomenologically as 

Ye{r) = a + 6arctan[(r - ro)/rs] , (6) 

with a = 0.24 and b = 0.165. The parameters ro and 
rg describe where the rise takes place and how steep it 



is, respectively. In the bottom panel of Fig. [T] we show 
two Ye(r) profiles for two different choices of these pa- 
rameters. The radius where A(r)Ye(^) crosses the hori- 
zontal bands wh ((^l) determines the well known H (L) 
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Mikheycv-Smirnov-Wolfcnstein (MSW) resonances 
Is^. For the A and Ye profiles shown in Fig. [1] and en- 
ergies typical in SNe the position of both resonances ru 
and tl he above 10^ km. 

The other non-zero element in Vgtd arises from ra- 
diative corrections to neutral-current z/^ and Vt scat- 
tering. Although there are no fi nor r leptons in nor- 
mal matter, they appear as virtual states causing a shift 
AV^r — V^GpY^^UB between i/^ and i^r due to the 
difference in their masses. It has the same effect on neu- 
trino dispersion as real r leptons with an effective abun- 
dance 
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= 2.7 X IQ-^ . 



(7) 

where rig = Up and Yn ~ 0.5 was assumed. Analo- 
gously to the H- and L-rcsonances, the radius where 
X{r)Y^^ sa wh or, equivalently A(7') ~ = ^a/Y.^^ , 
defines the /ir-rcsonance [g^l- The upper green band in 
the top panel of Fig. [1] represents uj^t- for energies be- 
tween 5 and 50 MeV. In the present analysis we have 
not considered the flavor non-universal radiative correc- 
tion for neutrino refraction in the presence of a neutrino 
background [67 1. 

In order to construct the term in the Hamiltonian de- 
scribing the non-standard neutrino interactions with a 
fcrmion / wc parametrize them with the effective low- 
energy four-fermion operator: 

Cnsi = ~ei''p2V2GFii^al^Li^p)ifrPf) , (8) 

where P = L, R and / is a first generation fermion: 
e, u, d. We neglect the possible NSI effect arising in 
J/ — J/ interactions [fi^. The coefficients e;^^ denote the 
strength of the NSI between the neutrinos of flavors a 
and /3 and the P— handed component of the fermion /. 
Throughout the article we assume neutrinos propagating 
across an unpolarized medium [g^. Therefore their evo- 
kition in matter will be affected by the vector coupling 
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constant of the NSI, e; 
J 



fv 

a/3 



-afS ' 



-ei^^,seeRei. [eg]. We 



also consider e^^ S 3?, neglecting possible CP violation 
in the new interactions''. Under these assumptions Vnsi 
can be expressed as, 



(Vnsi)a/: 



E 

f—e,u,d 



(9) 



where Nf represents the fermion / number density. For 
definiteness and motivated by actual models, for exam- 
ple, those with broken R parity supersymmctry we take 
for / the down-type quark. Therefore the NSI potential 
can be expressed as follows. 



(10) 



From now on we will not explicitly write the superindex 
d. In order to further simplify the problem one can always 
redefine the diagonal NSI parameters so that e^^ = 0, 
since subtracting a matrix proportional to the identity 
leaves the physics involved in the neutrino oscillation un- 
affected. Thus, in principle at least five of the six inde- 
pendent Eap parameters should be taken into account. 
Nevertheless, the exhaustive description of NSI in SN 
neutrinos developed in Ref. 36[ shows that the physics 



involved can be described in terms of e'^^ and e'^^ , which 
are just a suitable combination of e's. This motivates us 
to illustrate the interplay that could arise between col- 
lective effects and NSI by only considering two non-zero 
NSI parameters: Eer and Ett-, describing flavor-changing 
(FC) processes and non- universality (NU), respectively. 
Therefore the term in the Hamiltonian responsible for 
the interactions with matter can be written as 



/ 



V = A(r)(2-ye) 







V £e 






Er 







(11) 



The range of values for the NSI parameters we consider 
is for the off-diagonal term 10~^ < |£er| ^ few x 10"^. 
This prevents any significant NSI-induced reduction of 
the electron fraction Ye during the core collapse [tiI, [t^ . 
For the diagonal term we assume \£rT\ S, 0.1, allowed by 



the current experimental constraints 73l476l| 



Finally, the third term in the Hamiltonian accounts for 
the collective flavor transformations induced by neutrino- 
neutrino interaction. In spherical symmetry the EOMs 
can be expressed as a closed set of differential equations 
along the radial direction [45,|49|. We solve them numer- 
ically as described in Ref. [45|. The factor (1 — Vq • Vp) in 
the Hamiltonian, Vp being the velocity, implies "multi- 
angle effects" for neutrinos moving on different trajecto- 
We consider, though, the single-angle 



381 



41 



ries 

approximation by launching all neutrinos with 45° rela- 
tive to the radial direction. This approximation is valid 
for the neutrino fluxes assumed in this analysis provided 
that the neutrino density exceeds the electron density. 
The strength of the neutrino-neutrino interaction can be 
parametrized by 



Mo 



(12) 



where the fluxes arc taken at the radius R. Following 
Ref. 45| we shall assume fiQ — 7 x 10'^ km^^. In the 
single-angle approximation the radial dependence of the 
neutrino-neutrino interaction strength can be explicitly 
written as 



Mo 



2r4 



(13) 



In the top panel of Fig. [T] we show the typical /i(r) 
profile for yUg = 7 x 10^ km^^. One final property of 
SN neutrinos with important consequences for our study 
is the hierarchy of fluxes obtained in SN models. The 
typical conditions of the proto-neutron star lead to the 
following hierarchy of fluxes > F^^ > F^ , with 
I'x standing for i/^, i/,-, i/^ and Dr- As in Ref. 
we express the Icpton asymmetry with the parameter 
e = (F^l - F^^)/{F?^ - Fp,)- Throughout the analysis 
we shall consider e = 0.25. The equal parts of the fluxes 
drop out of the EOMs, so as initial condition we use in the 



monoenergetic case F^ 



Q&ndF^ ^ (l + e)Fp^. 



III. NON-STANDARD INTERACTIONS IN THE 
PRESENCE OF COLLECTIVE EFFECTS 



^ For the sake of simplicity we will omit the superindex V . 

* Possible CP effects in the case of SUSY radiative corrections on 

/i— r neutrino refraction including R-parity breaking interactions 

have been considered in Ref. [toI 



In the absence of NSI and collective effects the neu- 
trino propagation through the SN envelope is basically 
determined at the well-known MSW resonances defined 
in Secini They arise when the matter potential becomes 
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of the same order as the kinetic terms in the Hamiltonian. 
The L-resonance occurs always for neutrinos whereas the 
H-one takes place for (anti)neutrinos for (inverted) nor- 
mal mass hierarchy. For our matter profiles and the val- 
ues of t?i2 and T^ia both resonances are adiabatic 78, 7^. 
Moreover both involve electron neutrino flavor and hap- 
pen in the outer layers of the SN envelope, see top panel 

of Fig. m 

In addition, the /xr-resonance is also adiabatic 66| , but 



occurs between the and v^- or i/^ and Dr depending on 
the neutrino mass hierarchy and the 1923 octant. How- 
ever, when considering the neutrino self-interaction this 
resonance can also cause significant modifications of the 
overall z/g and De survival probabilities |53l | . According to 
the discussion in the previous section the /ir-resonance 
occurs at 



R 



R 



(14) 



Due to the smallness of the /ir-resonance happens at 
deeper layers than the H- and L-resonanccs. In particu- 
lar, for cjh = 0.3 km^^, = 0.5, and Aq = 5 x 10^ km^^ 
(4 X 10^ km~^) r^r = 770 km (71 km), see the intersec- 
tion between the upper green band and the profiles Ai(r) 
and A2(r) in the top panel of Fig. [TJ 

The consequence of the addition of an NSI term such 
as that of Eq. (|T0)) is twofold. First, it will affect the 
MSW resonances. For the values assumed here the main 
effect on the H- and L-resonances will be just a slight 
shift in the resonance point 36|, |80j ■ The consequences 



for the /iT-rcsonance can be more drastic. For sufficiently 
large values of \Stt\ a negative sign can change the res- 
onance channel, from to P or viccversa, depending on 
the octant of i?23- On the other hand, it can significantly 
modify the position of the resonance. In the presence of 
NSI the jLtr-resonance happens at 

1/3 



' fJ,T 



R 



Aon" 



off 
nsi 



1/3 



Ao 

UJ 



where we have defined 



-y^cff _ -y^cff 
r,nsi — T 



(2 - Ye)Srr : 



(16) 
(17) 



In particular, for {errl > Y^ /(2 — Ye) the value of wJJ^' 
decreases, and therefore r^,- is pushed outwards with re- 
spect to the standard case. The lower green band in the 



top panel of Fig. [T] represents lu^^ for Err = 0.1, Ye = 0.5 
and typical SN neutrino energies. For E = 20 MeV 
and the matter profile corresponding to A2.0, the po- 
sition of the /ir-resonance moves out to a radius of 
r^r w 1.3 X 10^ km. 

The second important consequence is that the new NSI 
terms can induce additional resonances [l7| . If we con- 
sider the inner layers, defined as those where r <^ th, 
in the absence of neutrino self-interaction, one may have 

« V, as given in Eq. (fTT|) . One can see that a novel 
resonance, which wc call /-resonance, / standing for "in- 
ternal" , will arise when the condition Hee = Hrr is sat- 
isfied. This occurs when the value of the is of the 



same order as the electron fraction Y^ [17[ . Current con- 
straints on the Eq^'s imply that small values of Y^ are 
required for these NSI-induced internal resonances to oc- 
cur. Hence this condition is only fulfilled in the most 
deleptonized inner layers, close to the neutrinosphere, 
where Ye reaches values of a few % 
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36|. If 



we neglect the contribution from Y°^ the corresponding 
resonance condition can be written as 

2£"rr 



YJ 



1 



(18) 



In the bottom panel of Fig. [T] we show as a horizontal 
band the range of Y/ required for the /-resonance to take 
place for Err < 0.1. One sees how for the Ye(r) profiles 
found in numerical simulations the resonance condition 
can only be satisfied in the inner layers. For typical values 
of Ye one expects to have the /-resonance for Err ^ 10^^. 
Moreover, as seen in [36[, the range of jeerl considered 
ensures adiabaticity. It must be noted that, in contrast to 
the standard H- and L-resonances, related to the kinetic 
term, neither the density nor the energy enter explicitly 
into the resonance condition, which is determined only 
by the electron fraction Ye. Moreover, in contrast to 
the standard resonances, the /-resonance occurs for both 
neutrinos and antineutrinos simultaneously 17 1. 



At the same time, also in the internal region, the neu- 
trino fiux emerging from the supernova core is so dense 
that, neutrino- neutrino refraction can cause nonlinear 
fiavor oscillation phenomena 
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The crucial effect 



is a collective mode of pair transformations of the form 
I'ef'e -T- Vxi^x, where x represents some suitable superpo- 
sition of ly^ and ■ This pair- wise form of flavor trans- 
formation leaves the net flavor-lepton number flux un- 
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changed. For the hierarchy of neutrino fluxes assumed 
this conversion occurs only when the neutrino mass hi- 
erarchy is inverted. Collective flavor transformations 
start after the synchronization phase, where ^(rgyn) ~ 
2a;H/(l — VT~Fe)^, and extends a few hundred km in 
the so-called bipolar regime until /i(rbip) ~ 42|. At 
larger radii /i(r) < and the neutrino self- interaction 
becomes negligible. For our chosen /xg, an excess flux 
of 25%, and wh = 0.3 km~^, we find a synchronization 
and bipolar radius of rsyn — 100 km and rbip — 330 km, 
as indicated in Fig. [1] by dark and light vertical gray 
bands, respectively. One important consequence of this 
flavor transformation in the context of three neutrino fla- 
vors is its potential sensitivity to deviations of i?23 from 
maximal mixing. As pointed out in Ref. |53|, in the par- 
ticular case that the /ir-resonance takes place outside 
the bipolar radius, r^r > T^bip, the final states Vxi^x, and 
therefore the final and survival probability, depend 
crucially on the octant of t?23- According to Eq. ([T5]) and 



the definition of rbip this condition is satisfied when 



> 



?'bip 

R 



ll/4 



~ LlxlO** km"' 
(19) 

In the standard case the possibility to discern the ??23 
octant only occurs for large density profiles, Aq > 7 x 
10^ km"^, i.e. at early times. This situation would corre- 
spond to the Ai (r) profile in top panel of Fig.dJ but not to 
A2 (r) . However the presence of NSI terms in the Hamilto- 
nian may shift the /ir-resonance to outer layers, making 
this condition more flexible. For instance, for Ett = 0.1 
the previous condition requires only Aq > 1.3 x 10^ km~^, 
see lower green band in the top panel of Fig. [TJ Therefore 
the presence of NSI could keep the possibility to distin- 
guish between the two 1923 octants for several seconds. 

The self-induced flavor transformations however do not 
occur for arbitrarily large density proflles. If the electron 
density Ue signiflcantly exceeds the neutrino density rin 
in the conversion region they can be suppressed by mat- 
ter 8l| . This is a consequence of neutrinos traveling on 
different trajectories when streaming from a source that 
is not point-like. This multi-angle matter effect can be 
neglected if in the collective region, prior to the synchro- 
nization radius, we have 



The limiting condition can be determined by imposing 
Eq. (pn|) at Tsyn. Taking into account Eqs. ^ and 
we obtain 



^'eC'^syiOAo < /iO 



(21) 



A'^(r)^ye(r)A(r)^«MO- 



(20) 



If we assume Yf, = 0.5 and Vsyn = 100 km then this con- 
dition amounts to Aq <C 1.4 x 10^ km"^. In the examples 
considered in the top panel of Fig. [1] this is fulflUed for 
A2,o- The condition \%{r) ^ /i(r) is then satisfled in the 
bipolar region. As a consequence collective effects are 
not matter suppressed. This is not the case of Ai(r). 
In the standard case the limiting Aq = 1.4 x 10^ km"^, 
above which multi-angle matter effects suppress the col- 
lective effects, is, though, smaller than the minimum 
Aq = 7 X 10^ km"^ reqiured for the effect to be im- 
portant. The situation could drastically change in the 
presence of NSI. As previously discussed, non-zero NSI 
non-universal (diagonal) parameters could help moving 
the /tr resonance out of the rgyn even for Aq smaller than 
1.4 X 10^ km"^. The consequence is that large enough 
[e^T-l would make the neutrino propagation through the 
SN envelope highly sensitive to the 7923 octant. 



IV. CLASSIFICATION OF REGIMES 

In this section we summarize all the information for- 
merly introduced. Taking into account the conditions 
given in Eqs. ([TH]), ([11]), and (HJ) we can roughly 

identify four different regimes of the neutrino propaga- 
tion depending on Aq and Et-t- This scheme is displayed 
in Fig. H for e = 0.25 and /io = 7 x 10^ km"!. It is 
important to notice that the size of the different regions 
depends on the particular values of 2e and /zq- Never- 
theless, unless kinematic dccoherence between different 
angular modes is triggered [4^, this schematic picture 
certainly holds. 

To first approximation the four regions can be defined 
in terms of matter suppression (or not) of collective ef- 
fects and presence (or not) of the internal /-resonance. 
Equation ((2T|) is depicted as a horizontal solid line at 
Aq = 1.4x 10^ km"^. For higher Aq matter suppresses col- 
lective effects whereas for smaller densities collective ef- 
fects are present. For intermediate values, \*{r) A'(r), 
there would be a matter induced decoherence 



81|- To 
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make the discussion as simple as possible we will only 
consider the extreme cases. On the other hand the ver- 
tical dashed line at Err — 10"'^ indicates the presence 
(right) or absence (left) of the NSI-induced /-resonance. 
This simple scheme becomes further complicated if one 



no I - resonance 




FIG. 2: Different regimes of the neutrino propagation depend- 
ing on the value of Ao and e-rr, as described in the text. 



adds the possibility that the //r-resonance lies outside 
the bipolar region. In the next subsections we analyze in 
detail the different possibilities. 

A. Region I 

On the upper left corner we have the region I, defined 
by Aq > 1.4 X 10"^ km^i and e^^ < 10"^. Accord- 
ing to the previous discussion, this range of parameters 
leads to no collective effects, since they are suppressed 
by matter, and no /-resonance. Assuming that the L- 
and /f-resonances are adiabatic the v^, and survival 
probability is then only fixed by the mass hierarchy. The 
NSI terms will lead at most to a small shift in its posi- 
tion 
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the presence of a dip in the uppermost curves at high den- 
sities, corresponding to the T4e entry in Eq. (jlip. This is 
a direct consequence of the radial dependence of Yeif), 
as given in Eq. ^ . The arrows in these figures represent 
the transitions caused by the collective effects, and must 
therefore be ignored when these are not present. In the 
normal hierarchy case v^, and leave the SN as v-^ and 
whereas for inverted mass hierarchy they escape as 
vi and vj, for any octant. The survival probabilities can 
then be written as P(ye -> « sin^ iJis (sin^ "dyi) and 
Pipe i^e) = cos^ ??i2 (sin'^ T^is) for normal (inverted) 
mass hierarchy. Figure [4] represents in solid lines the ra- 
dial evolution of Pee a-nd assuming Ag = 1x10^ km~^, 
WH = 0.3 km-\ sin^ d-ij. = 0.5, and e = 0.25. The verti- 
cal bands indicate where the resonance conversions take 
place. In order to perform the plot we have artificially 
set po = 0. We want to remind here that both and 
Pee are normalized to the flux, and therefore, while pee 
corresponds directly to survival probability, pee must 
be corrected by a factor (1 + e) in order to obtain the cor- 
responding survival probability, pee = P(ye — > i^e)(l + e)- 



B. Region II 

The region II, on the upper right corner, is defined 
by Ao > 1.4 X 10'^ km^i and e^r > 10"^. As in re- 
gion I the matter density is so high that prevents neu- 
trinos from undergoing collective effects. However, the 
values of the diagonal NSI terms in this region are large 
enough to fulfill Eq. p8|) , causing the /-resonance to ap- 
pear. In contrast to the previous case and Ve are now 
created as v^^ (y™) and (yf ) for normal (inverted) 
mass hierarchy ^ . They cross adiabatically all resonances 
and leave the SN as {v\) and Vi (y\) for normal (in- 
verted) mass hierarchy [36[ . In Fig. [5] we show the level 
crossing scheme analogous to Fig. [3] but in the presence 
of the /-resonance. The survival probabilities are now 



In Fig. |3]we show the level crossing scheme for normal 
(top), and inverted mass hierarchy for 1923 < 7r/4 (mid- 
dle), and z?23 > 7r/4 (bottom). It is remarkable to note 



The different neutrino basis used in tfie analysis are: Va witfi a = 
e, /i, T representing the flavor basis; Vi with i = 1, 2, 3 standing 
for the mass basis, diagonalizing the Hamiltonian in vacuum, 
_ff = f2; the matter basis, v™- , diagonalizing the Hamiltonian in 
absence of neutrino self-interactions = f2 -|- V; and the basis 
v'^ defined as the one diagonalizing the Hamiltonian in vacuum 
in the case that all mixing angles are very small {after rotating 
the matter term away). 
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FIG. 4: Radial dependence of density matrix elements pee and 
Pee corresponding to regions I (solid) and II (dashed) in Fig. [2] 
Top panel represents normal mass hierarchy and bottom panel 
inverted mass hierarchy. We assume 

ujn = 0.3 km~^, and sin^ 1923 = 0.5. Vertical bands indicate 
regions where resonances take place. 



FIG. 3: Level crossings in the absence of /-resonance for 
normal and i923 < 7i"/4 (top), and inverted mass hierarchy 
for i923 < 7r/4 (middle), and i?23 > 7r/4 (bottom). The 
dashed and solid arrows in the middle and bottom indicate 
the pair transformations due to collective effects happening 
after (dashed) or before (solid) the pr-resonance. 



P{Uf. V(.) = P(9e ^e) ~ sin^ 'i?i2 (cos^'i?i2) for nor- 
mal (inverted) mass hierarchy. The black dashed lines in 
Fig. |4] show the expected radial evolution of pee and pee, 
respectively, when neutrinos and antineutrinos undergo 
an adiabatic /-resonance. The width of the green band 



accounts for the oscillatory behavior of neutrinos. As for 
region I, we have made the calculation assuming /ig = 0. 
However, we have analyzed the single energy and multi- 
angle case within two-flavor framework for the range of 
parameters here discussed, and have verified that collec- 
tive effects are indeed suppressed and the /-resonance is 
present for both normal and inverted hierarchies. This 
means that the behavior in region II corresponds indeed 
to the case discussed in Ref. 361. 
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FIG. 5: Same as Fig.[3]in the presence of /-resonance. Dashed 
arrows represent pair transformations due to collective effects 
occurring after the /-resonance. In the bottom panel there are 
two possibilities for this transition, corresponding to an adi- 
abatic or non-adiabatic /-resonance. The solid arrows in the 
middle and bottom panel refer to collective effects happening 
before the /-resonance. 



C. Region III 

Let us now consider the lower part of Fig. [2l i.e. when 
'^o ^ 1-4 X 10^ km^^. The main feature of this sce- 
nario is the presence of collective effects. As discussed in 



Ref. 53[ , and here reviewed, these in turn depend on the 
relative position of the /xr-resonance with respect to the 
synchronization and bipolar radius. We can then distin- 
guish two different regimes: On the bottom left corner 
we define region III by the condition r^.r ^bip, and on 
the bottom right corner we have region IV defined by 
r^T- > J'bip. In the middle of both there is a transition 
region whose width we have determined numerically. 

Let us now discuss region III. According to Eq. 
this range of parameters satisfies the condition 



\Q\Yf + (2 - Y,)err\ < 1.1 X 10* km" 



(22) 



which, for Ye = 0.5, roughly amounts to \q£tt 
7.3 X 10'^ km^^, see Fig. [21 This situation can be re- 
duced to the standard two-flavor scenario. In order to 
better understand the consequences of collective effects 
it is convenient to realize that the impact of ordinary 
matter can be transformed away by going into a rotat 
ing reference frame for the polarization vectors 
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Collective conversions proceed in the same way as they 
would in the absence of the matter term T4td, except 
that the effective mixing angle is reduced. After this 
rotation the mass eigenstates v[ now approximately co- 
incide with the interaction eigenstates v^. In particular, 
the initial states and can therefore be essentially 
identified with i/( and v'^, respectively. If the neutrino 
mass hierarchy is normal, we begin in the lowest-lying 
state and no collective effects take place, see bottom left 
panel of Fig. [S] (also in Fig. 4 of Ref. (Sg). The situ- 
ation is then similar to that in region I, see top panel 
of Fig. [31 However in the case of inverted mass hierar- 
chy both v'-y^ and v'^ correspond to the intermediate state. 
The effect of the self-interaction is to drive them to the 
lowest-lying states, which in this case are v'^ and v'^^ see 
top left panel of Fig.[6l In terms of matter eigenstates vf" 
this is shown by dashed arrows in the middle and bot- 
tom panels of Fig. [3l In the case of a fraction equal to 
tFy^ is not transformed and stays in v^^ and evolves as 
in the absence of neutrino-neutrino interactions, i.e. adi- 
abatically through the L-resonance. The rest of z^e are 
transformed to . As a result the final v^, flux, normal- 
ized to the initial v,, flux, is expected to be approximately 
pflnai ^ esm^'dx2 + sin^t^ia ~ 0.08. On the other hand, 
after the pair transformation v^. cross the i?-resonance 
adiabatically and leave the star as leading to a final 
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FIG. 6: Neutrino vacuum level diagrams without NSl effects 
(left panels) and with the NSI-induced I resonance occurring 
before the bipolar conversion (right panels). The top (bot- 
tom) panels correspond to inverted (normal) neutrino mass 
hierarchy. In each panel the two possible i?23 are indicated. 
The 12 and 13-mixing angles are assumed to be very small, 
mimicking the effect of ordinary matter. The effect of collec- 
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normalized flux of approximately p^" ^' = cos^ i?i2 — 0.7. 
This can be seen in Fig. [71 where we show in solid lines 
the radial evolution of ly^ and Dg for inverted mass hier- 
archy assuming Aq = 4 x 10^ km^^, cjh = 0.3 km^^, and 
sin^ i?23 = 0.4 (top) and 0.6 (bottom). As can be seen in 
the figure, the result is independent of the ^23 octant. 



FIG. 7: Radial dependence of pee and pee for region III with 
Ett = (solid) and IVa (dashed) with Ett = 3 x 10"'^ for 
inverted mass hierarchy, and sin^ i?23 = 0.4 (0.6) in the top 
(bottom panel). In both cases Ao = 4x 10^ km~^ and Ser = 0. 
The bands around the lines represent modulations. Vertical 
gray bands stand for synchronized (dark) and bipolar (light) 
regime. Resonance regions are also displayed. 



D. Region IV 

Finally, neutrinos with parameters in the right bot- 
tom corner (region IV) will feel both collective and NSl 
effects. This region of parameters is defined by the condi- 
tion that the /zr-resonance lies outside the bipolar region. 
According to Eq. (fTg]) this amounts to 

Ao[r;^ + (2 - Ye)err] > 1-1 X 10^ km~^ . (23) 

As discussed above, for the standard case and = 0.5 
this is satisfied for Aq > 7 x 10® km^^, which implies 
a strong matter suppression of the collective effects, see 
Fig. 121 However, if NSl diagonal parameters are of the 



order of \err \ > 7.3 x IOVAq (km"^) then one can avoid 
the matter suppression condition. Therefore the first NSl 
effect is to increase the value of Y^^^^ so that the Aq re- 
quired to have the /zr-resonance outside rbip is still con- 
sistent with the presence of collective effects. Moreover, 
if Srr is of the order of a few % the condition given in 
Eq. ([18)) is fulfilled for the typical values of found in 
SNe. Thus, in the region IV we can distinguish two sub- 
sets of parameters denoted by IVa and IVb defined by 
the absence or presence of the /-resonance, respectively. 
For defineteness we set the boundary at Et-t = 10~^. 

Let us first consider the IVa region. Depending on Aq, 
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i.e. on the instant considered, this range of parameters 



imphes values of 



from 10 



to 



10~2. Ahhough 



values are not high enough to induce the / resonance they 
are sufficiently large to push the /ir-resonance outside the 
bipolar region. The situation is therefore analogous to 
the one described in Ref. [5^. That means a flavor pair 
transformation Vf.^!, — >■ v^v^ due to collective effects only 
for inverted neutrino mass hierarchy, like in region III. 
However, the final matter eigenstates depend on the 1^23 
octant. In the middle panel of Fig. |3] we show with solid 
lines the pair conversion for 1^23 in the first octant. In 
terms of matter eigenstates, Ve and i^e are transformed 
into I'™ and respectively. The presence of the ^r- 
resonance in the neutrino channel leads to a difference of 
with respect to region III. In the top panel of Fig. [7] we 
show with dashed lines the evolution of as function of 
the distance, for 10 ^ and sin^ i?23 = 0.4. In the 

collective bipolar conversions, the excess e of Ve over Vf, 
remains as vl^ whereas the rest will be transformed to v^. 
As a consequence, the original z/g flux leaving the star can 
be written as pf"''' = esin^ ??i2 + cos^ i9i2, which in our 
particular case amounts to roughly 0.75. If 1^23 belongs to 
the second octant the /xr-resonance takes place in the an- 
tineutrino channel, see bottom panel of Fig. [31 The pair 
Ue and i^e is driven to the lowest-lying states, which in 
this case are and for neutrinos and antineutrinos, 
respectively. Therefore, for Vf, the situation is completely 
analogous to that in region III, whereas Ve leave the star 
as V2. The radial evolution of for Sfp — - 3 X 
and sin^ 1^23 = 0.6 is displayed with dashed lines in the 
bottom panel of Fig. [71 It is remarkable that neither i/^ 
nor Pe undergo the H resonance, and therefore are blind 
to the possible effect of the outwards propagating shock 
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wave 

It is important to notice that the same effect observed 
for the different octants of ^23 can be obtained by fixing 
the octant and changing the sign of Err- This can be 
easily understood if we study the /ir-resonance condition, 

A(r) [Yf + (2 - Ye)err] ^ -UJR cos2 cos 2Z?23 , (24) 

where we have neglected subleading solar terms. This 
condition dictates the channel where the resonance takes 
place. In the standard case this is determined only by 
the hierarchy of neutrino masses and the octant of '!?23- 
In the presence of NSI, though, the sign of left-hand side 



of the equation depends on that of Ett, what therefore 
affects directly the resonance condition. As a result, the 
same result of Fig. [71 is obtained by changing the sign 
of Err and the octant of 'i?23, i-c Err = — 3 X lO^'' and 
second octant for the panel on top and first octant in 
the bottom panel. In summary, in the presence of non- 
universal NSI parameters with {EttI ^ 10"'^, SN neutri- 
nos are sensitive to the octant of i?23, the absolute value 
of the non-universal NSI parameter as well as its sign. 

Finally, for higher values of the non-universal NSI 
parameters, Err ?J 10^^ (region IVb), the internal I- 
resonance will arise. In this case one must analyze the 
interplay between collective effects and the /-resonance. 
This is discussed in the next section. 



V. COLLECTIVE EFFECTS AND 
NSI-INDUCED /-RESONANCE 



In this section we analyze region IVb, defined by Aq 



< 



1.4 X 10^ km^^ and Err ^ 10^^, where both an adiabatic 
/-resonance and collective effects are present. If the /- 
resonance is not adiabatic then neutrinos within region 
IVb evolve exactly as in IVa. 

As can be inferred from Fig.[Tl one of the main features 
of this scenario is that both effects happen nearly in the 
same region, namely the deepest layers right above the 
neutrinosphere. This means that the final result will also 
depend on the relative position between the bipolar re- 
gion and the location of the /-resonance. Schematically 
two extreme scenarios can be identified. In one case the 
rise in the Ye, and consequently the /-resonance, takes 
place before the bipolar conversion region, see Y^ in the 
bottom panel of Fig. [H In the second scenario one has 
first the bipolar conversion and then neutrinos traverse 
the /-resonance, see Y^ in the bottom panel of Fig. [H 

A. First NSI /-resonance 

Let us analyze here the particular case where the /- 
resonance happens in deeper layers than collective effects. 
This situation corresponds to the Y^ in the bottom panel 
of Fig. [H and, according to SN numerical simulations it 
is the most likely situation. 

The main consequence of an /-resonance in the inner- 
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most layers, right after the neutrinosphere, is an inversion 
of the neutrino fluxes entering the bipolar region. For the 
initial flux pattern assumed in Sec.|lT]this implies the fol- 
lowing new pattern after the /-resonance: Fi,^ = = 0, 
F,^^ = 1 + e, and Fp^ = 1 normalized to Ff^^ . In contrast 
to the standard case, under this condition collective ef- 
fects arise in the case of normal mass hierarchy. 

This can be understood using the pendulum anal- 
ogy [42l | in the corresponding reduced two flavor scenario. 



and keeping in mind that the bipolar conversion drives 
the neutrinos to the lowest-lying states. In the normal hi- 
erarchy the system is already created near the minimum 
of the potential. Thus, in absence of the /-resonance col- 
lective effects are not present. However once an adiabatic 
/-resonance occurs the neutrino flavor is swapped and the 
system is driven to the maximum of the potential. The 
original Ve (i^e) are transformed to i^r (i^r) and, in terms 
of the vacuum level diagrams, occupy now the highest 
states right before collective effects are switched on. In 
the right panels of Fig. |6] we show the vacuum level di- 
agrams after an adiabatic /—resonance has occurred. In 
the case of normal mass hierarchy (bottom right panel) 
bipolar effects act taking 1/^ (I'e) to the lowest-lying states 
j/J (f'i)i which in this case do not depend on the 7?23 oc- 
tant. This bipolar conversion can be also seen in term of 
matter eigenstates as dashed arrows in the top panel of 
Fig. [51 As a result both the /-resonance and the induced 
collective effects basically cancel each other. In the top 
panel of Fig. |8] we show the radial evolution of i^e and 
De. This cancellation between the /-resonance and col- 
lective effects is complete for z/g, which leave as Di, but 
not for z^e! its excess e over D,, remains as 1/3 (1^2) for '^23 
in the first (second) octant. Therefore for a monochro- 
matic neutrino flux we obtain, after propagating through 
the outer resonances, p^^'^^ = esin^ i9i2 + sin^ 1^13 « 0.08 
for both 1^23 octants, see top panel of Fig. |8l instead of 
simply sin^ i?i3 = 10^^ as in top panel of Fig. |4l Hence, 
by comparing the top panels of Figs. |4] and [8] one real- 
izes that, except for the excess e, the situation for normal 
mass hierarchy is basically the same as in regions I, HI, 
and IVa. 

While this is true in the monoenergetic case a specific 
signature can be observed if we do not restrict ourselves 
to that case but consider the whole energy spectrum. 
The top panel of Fig. [HI displays the and z/^ fluxes at 



the neutrinosphere, and . We have assumed the 
parameterization given in Ref. 



e{E) = c,^ 



E 



exp 



E 



(25) 



with {EpJ ^ 12MeV, {E^J = 15MeV, {E^J = IBMeV, 
= 5, pjy^ ~ 4.5 and /?^^ = 4. The normalization 
Cj/^ has been chosen such that F^ = J ff {E)dE = 1, 
/^R = 1 -t- Ke and F^^^ = 1 - k, with k = 0.15. 

In the bottom panel of Fig. [HI we show the v^. fluxes 
after the bipolar region. In solid dark red lines we repre- 
sent the case of normal mass hierarchy in region IVb. By 
comparing the two panels one sees how the conversion 
^ i^x takes place only at low energies. This is exactly 
the opposite of what happens for the standard case (in- 
verted mass hierarchy in region HI) , shown as solid light 
red lines, where the untransformed flux concentrates at 
low energies 46|, |47[ . For completeness we show also the 



other cases. In the region of parameters I, HI and IVa 
there is neither collective effects nor /-resonance for nor- 
mal hierarchy, then the fluxes after the bipolar region 
coincide with the initial ones, f,j^ = f^^ and /j^^ = . 
In region II, the /-resonance implies a complete conver- 
sion — > which leads to a spectral swap, f,j^ = 
and /^^ = 

The case of inverted mass hierarchy is more subtle. Ac- 
cording to the previous discussion one would expect no 
collective effects after neutrinos traverse the /-resonance. 
The system starts its evolution near the maximum of the 
potential and, in the absence of NSI, the bipolar con- 
versions would take it to the minimum. What the /- 
resonance is doing in this language by swapping the flavor 
eigenstates is to take the system to the minimum of the 
potential before any collective effects can arise. The new 
stable situation prevents bipolar conversions, leaving the 
system unchanged until the outer resonances are reached. 
And this is indeed what happens if i?23 lies in the first 
octant. In terms of the vacuum levels diagram the origi- 
nal Ve (i^e) are transformed after the /-resonance into 
{Pt), which are already the lowest-lying states, see the 
top right panel in Fig. [6l Hence no collective effects take 
place, and, as can be seen in the middle panel of Fig. [5] 
after traversing the outer resonances and leave the 
star as and The radial evolution is shown with 
solid lines in the bottom panel of Fig. [51 
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FIG. 8: Same as Fig.|4]for region IVb, for normal (top panel) 
and inverted (bottom panel) mass hierarchy. In the top panel 
it is shown the case sin^ i923 ~ 0.4, while in the bottom panel 
the case sin^ ^23 = 0.6 is also displayed. In both cases Ao ~ 
4 X lO** km~^ and e^r = 5 x 10"^ and e^^ = 10"^. The 
/-resonance is assumed to occur before the collective effects 



However if 'i?23 lies in the second octant things are dif- 
ferent. As in the previous case, the /-resonance converts 
neutrinos created initially as and De into Vj- and Dr, 
respectively. But these states are now the highest states, 
1^2 and This means, that in contrast to the first-octant 
case, when the neutrinos traverse the bipolar regime they 
will be driven to the lowest-lying states, i.e. 1/21^2 — >■ 1^31^3, 
see top right panel of Fig. [6l In the bottom panel of 
Fig. [8] we show with dashed lines the evolution of ly^ and 
De as function of distance for sin'^ i?23 = 0.6. The bipolar 
conversion can not be directly seen as it occurs between 
states containing explicitly neither v^. nor Dg. At the end, 
as shown in the bottom panel of Fig. [51 the original 
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FIG. 9: Top: and i/^ fluxes as emitted at the neutri- 
nosphere. Bottom: z^e flux after the bipolar region for differ- 
ent cases. In the region IVb it is assumed that the /-resonance 
occurs before than the bipolar conversion. Note a depletion 
in the flux of low-energy neutrinos, instead of usual spectral 
split. 



leave the SN as P2, see bottom panel of Fig. Hi In the case 
of the excess e over remains as 1^2 whereas the rest 
is transformed to V3. Therefore for a monoenergetic flux 
we find pfe'^' = ecos^ i?i2 -l-sin^ dis « 0.2. By comparing 
the bottom panel of Figs. |8| and [7| one realizes that this 
case is analogous to IVa. Since the collective effects do 
not affect i^e and Ve directly, considering neutrinos with 
an energy spectrum, one expects simply a complete swap 
of spectra, f^^ = /^^ and f^^ = /^^ after the bipolar 
region, like in scenario II, see bottom panel of Fig. |9| 

The final conclusion is that the propagation of neu- 
trinos with parameters in region IV is practically inde- 
pendent of whether the /-resonance is present (IVb) or 
not (IVa). The main consequence of the /-resonance is 
to swap the spectrum at low energies and only in the 
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case of normal mass hierarchy, in contrast to the trans- 
formation of high-energy v^^ happening for inverted mass 
hierarchy in the absence of /-resonance. 

B. First Collective 

For completeness we have also considered the possibil- 
ity that the bipolar conversion takes place before neutri- 
nos traverse the /-resonance. This situation corresponds 
schematically to the Y^^ profile in the bottom panel of 
Fig.ffl 

The case of normal mass hierarchy is completely analo- 
gous to the one of region II, that is, absence of collective 
effects and /-resonance. The i^e and are created as 
I/™ and P™, respectively. Therefore, if all resonances in- 
volved, /, /XT, //, and i, are adiabatic then they leave 
the SN as V2 and V2 , respectively, see top panel of Fig. [101 
The result is identical to the one shown with dashed lines 
in the top panel of Fig. |4l 

The situation with inverted mass hierarchy depends 
significantly on the "Qiz octant. Rotating the matter term 
away z^e and v^. are created as the intermediate states 
v'^ and v'y- Collective effects drive them to the lowest- 
lying states and v'^. However the corresponding mat- 
ter eigenstates are different depending on whether "diz 
belongs to the first or to the second octant. In the first 
case, most of Vf^ and Vf, end up as and before 
crossing the /-resonance, see solid arrows in the middle 
panel of Fig. [S] The excess e of v^, stays as vf . As a 
consequence, the final z^e and Vp^ fluxes, normalized to 
the initial i^e one, are = e cos^ i?i2 + sin^ t?i2 and 

Pee^' = sin^i?i3, respectively. See solid lines in the bot- 
tom panel of Fig. [TOl Except for the excess e in the 
net result is a cancellation of the collective effects and 
the /-resonance, leading to a similar result as in region I 
(solid lines in bottom panel of Fig. |4|). Qualitatively, the 
main difference shows up in the spectrum. The initial 
collective effects induces a "standard" swap only at high 
energies. Nevertheless, as neutrinos cross the /-resonance 
this split turns into an inverse one, with a swap at low 
energies. The final result right after the /-resonance is 
analogous to the case of normal mass hierarchy and the 
/-resonance happening first, discussed in Sec. IV Al and 
displayed with dark red solid lines in the bottom panel 
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FIG. 10: Same as Fig. [8] but assuming tliat the collective 
effects take place before the /-resonance. 

of Fig. [91 

If ')?23 lies in the second octant then most of Vf, and 
end up as i/™ and t?™. As can be seen in the solid 
lines in the bottom panel of Fig. [5l these neutrinos will 
not traverse the /-resonance, except the excess e of z/g, 
which stays as i^f. These neutrinos will be basically blind 
to the /-, //-, and L-resonances. The final fluxes will be 
therefore p^" = sin^ t9i3+e cos^ -0x2 and = sin^ ■dx2- 
This case is represented with dashed lines in the bottom 
panel of Fig. [TOl In the end the final evolution turns out 
to be similar to that in region IVb. 

VI. DISCUSSION 

In the previous sections wc have studied the conse- 
quences of NSI on the neutrino propagation through the 
SN envelope taking into account the presence of a neu- 
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trino background. We have analyzed the different situ- 
ations in terms of the non-universal NSI parameter Ett 
and the density at the neutrinosphere Aq. Depending on 
their values we were able to identify four extreme regions 
of the parameters where the evolution of the neutrinos 
have a specific pattern. 

In a realistic situation, though, we expect to find a 
combination of these situations depending on the instant 
considered. As mentioned in Sec.|lT]one expects the value 
of Ao to decrease with time as the explosion goes on. 
Therefore it is important to look for a time dependence 
in the neutrino propagation for given values of the NSI 
parameters. In Fig.[Tl]we show the composition in terms 
of mass eigenstates of neutrinos created as Ve (upper half 
of the panels) and (lower half of the panels) at the neu- 
trinosphere, when they leave the SN, as function of time 
for different neutrino mass and mixing schemes and for a 
given value of Err- The evolution in time shown in each 
panel is equivalent to consider Fig. [2l fix a value in the x 
axis corresponding to some Err and following vertically 
towards lower values of Aq. Depending on Ett and the 
instant considered one can distinguish different regions 
separated by vertical bands denoting transition phases. 
The position and the size of these transition bands are 
not constant but depend on time, as /iq and e do. Nev- 
ertheless, unless multiangle decoherence is triggered (by 
e.g. a strong reduction of e 01)' the sequence of differ- 
ent regimes undergone by the neutrinos is not expected 
to change drastically. 

Let us first discuss the antineutrino case. In the left 
panel we consider the standard case, i.e. Err = 0. On the 
left we have early times (or large Aq), which corresponds 
to the region I in Fig. [5] From the previous discussion, 
we know that for such a case and normal mass hierar- 
chy (upper box) an antineutrino created as De leaves the 
star as Di whereas in inverted mass hierarchy they es- 
cape as 1/3 due to the adiabatic iJ-resonance for both 
first and second ^23 octants (middle and lower boxes, re- 
spectively), see Figs. [3] and m At later times, Aq becomes 
smaller and matter can not suppress collective effects any 
longer, i.e. neutrinos enter region III. These affect only 
in the inverted mass hierarchy case "canceling" the H- 
resonance conversion and making the initial i?e escape as 
Pi. There is then a time dependence in the survival De 
probability for inverted mass hierarchy but not for the 



normal one. As can be seen in the panel this behavior 
does not depend on the z?23 octant, see Figs. |3] and [71 In 
terms of De survival probabilities there is then a transi- 
tion from sin^ 7?i3 w at early times to cos^ z?i2 ~ 0.7 
at later times, the details depending on the specific time 
evolution of A(r). 

Let us now take the middle panel with Err = 3 x 10~^. 
The situation at early times is the same as in the previous 
panel, described by the region I. However at intermedi- 
ate times the situation changes in the case of inverted 
mass hierarchy. Now the NSI parameters make the evo- 
lution go through region IVa before entering eventually 
region III. The /ir-resonance is pushed outside the bipo- 
lar region and then the degeneracy between the two ^23 
octants is broken: for '(?23 in the first octant leaves as 
Di whereas for the second octant they escape as i>2, see 
Figs. O and [T] At later times Aq further decreases and 
the ^r-resonance contracts to deeper layers within rgyn- 
Hence neutrinos cross to region III and the '!923 octant 
degeneracy is restored. Concerning the De survival prob- 
ability, as before, there is a transition from sin^ di3 w 
to cos^ 1)12 ~ 0.7. For 1923 in the first octant it is di- 
rect, whereas if 1923 fies in the second octant it goes from 
sin^ 'di3 « through sin^ ^12 ~ 0.3. As it was discussed 
in Sec. IIVDI an analogous effect would arise for a fixed 
■d23 and different signs of Err- For NH there is no time 
dependence. 

Finally, we consider the case where the NSI param- 
eters are large enough. Err = 5 x 10^^, to induce the 
/-resonance, right panel of Fig. 1111 Now, at early times 
neutrino propagation follows the prescription given in re- 
gion II. For normal mass hierarchy the original leave 
the star as D2, whereas for inverted they escape as Pi 
for both octants, see Figs. [S] and [5] After this phase 
neutrinos enter the region IVb and collective effects are 
switched on. For NH they cancel the /-induced resonant 
conversion, whereas for IH they break the degeneracy of 
the '^23 octant for inverted mass hierarchy. 

The bottom line is that if \£tt\ ^ a 10^'^ neutrinos 
cross the region IV during some seconds, and this could 
help disentangle the 1923 octant. If the octant were known 
one could obtain information about the sign of the non- 
universal NSI parameters. 

In the upper half part of the panels we show the same 
kind of plots but for neutrinos. The main difference with 
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respect to antineutrinos is that in the presence of col- 
lective effects Vf, are not fully converted like Vf^. Some 
fraction of them, corresponding to the excess over j>e , re- 
mains unaffected. This excess is represented in Fig.fTTjas 
a small portion at the right hand side of the correspond- 
ing boxes. As discussed, this excess of translates to a 
spectral split, so that the flavor spectral swap happens 
only for some energies. Whether these correspond to the 
low-energy tail or high-energy tail of the initial spectrum 
depends on the neutrino oscillation parameters, see bot- 
tom panel of Fig.[9l Therefore one could hope to use this 
additional information to break possible degeneracies be- 
tween different mass and mixing schemes and different 
values of the NSI parameters. It is important to note 
that the results here presented are based on the assump- 
tion of a particular choice of the initial neutrino fluxes, 
with the hierarchy > F^^ > F^^ . However, as shown 
in Ref. [86| , a different pattern for the initial fluxes could 
lead to a more complicated structre of the final spectrum 
with multilpe spectral splits. 

Throughout the analysis we have concentrated on the 
case of NSI with d quarks. Nevertheless most of the re- 
sults here presented can be generalized to the case of u 
quarks and e. In the first case the only effect is to shift 
the position of the /-resonance, since the resonance con- 
dition is modified to F/ — Srrl^ — £tt) [s^- ^or the 
case of NSI involving electrons the /-resonance is absent. 
But nevertheless its contribution to increase the value of 
y^^gj would also make the neutrino propagation highly 
sensitive to the i?23 octant and to its own sign, exactly 
as in the case of d quark. 

Last but not least we briefly comment on the possi- 
bility to observe the different regimes analyzed. This 
possibility will be hampered by several uncertainties in- 
herent in SN neutrinos. One is the lack of knowledge on 
the exact matter profile traversed by the outgoing neutri- 
nos. In our study we assumed a simple power law given 
in Eq. ([5]). This density profile will be significantly dis- 
torted by the passage of the shock wave responsible for 

. One of the 



in the spectra could further help disentangle between the 
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ejecting the whole SN envelope 
main effect will be to destroy the adiabaticity of the H- 
resonancc, which was assumed in our study. This effect, 
though, is not always present but depends on the neu- 
trino parameters. Therefore, far from being a problem, 
the time and energy dependence modulation introduced 



different scenarios here considered 

Another important source of uncertainties is our igno- 
rance of the exact initial fiuxes /^(/?). Although the 
initial fluxes during the first stage of the explosion, the 
neutronization burst, are rather model independent [oot . 
the expected number of events is very low. Most of the 
signal is generated later, during the accretion and cooling 
phases. The spectral features observed in the numerical 
simulations depend strongly on the properties of the SN. 
It is therefore necessary to set up strategies combining 
different observables to be able to pin down the underly- 
ing neutrino properties independently of the initial fluxes. 
For example spectral modulations are expected i f neutri - 
nos cross the Earth before being detected 
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One may also take advantage of the time dependence of 



the matter proflles A(r) and Yeir) themselves j36| . 



VII. SUMMARY 

Here we have considered the effect of non-standard 
neutrino interactions on the propagation of neutrinos 
through the SN envelope. We adopted a realistic three- 
neutrino framework, properly taking into account the 
presence of a neutrino background. We have found that 
for non- universal NSI parameters exceeding 10"^, i.e. 
[e^^l > 10^'^, the neutrino propagation becomes for some 
time sensitive to the "(^23 octant and the sign of Err- In 
particular, for > 10~^ an internal /-resonance may 
arise independently of the matter density. For typical 
values found in simulations this takes place in the same 
dense-neutrino region above the neutrinosphere where 
collective effects occur, in particular during the synchro- 
nization regime. This resonance may lead to an exchange 
of the neutrino fluxes entering the bipolar regime. The 
main implications are (i) bipolar conversion taking place 
for normal neutrino mass hierarchy and (ii) a transforma- 
tion of the flux of low-energy i/g, instead of usual spectral 
swap. 
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